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The relaxation of a dissipative system to its equilibrium state often shows a multiexponential 
pattern with relaxation rates, which are typically considered to be an intrinsic property of the 
system, independent of the initial condition and given by the spectrum of a Hermitian operator to 
which the initial Fokker-Planck operator is transformed by a similarity transformation. We show that 
this is not always the case. Considering the exactly solvable examples of standard and generalized 
Levy Ornstein-Uhlenbeck processes with a-stable initial distributions we show that the relaxation 
rates belong the spectrum of the corresponding quantum harmonic oscillator Hamiltonian only if the 
initial distribution belongs to the domain of attraction of the stable distribution defining the noise. 
Thus, in case of the standard Ornstein-Uhlenbeck process, broad a-stable initial distributions show 
a different relaxation pattern, and this pattern can persist as a long transient even for truncated 
stable initial distributions. 
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The relaxation of a physical system, prepared in a 
non-equilibrium state, to the equilibrium often shows 
a multiexponential pattern with decrements of single 
exponentials defining the relaxation rates. The same 
rates usually govern the decay of correlation functions 
of observables measured at different times in the 
equilibrium state. The rates are often considered as 
an intrinsic property of the system, independent of 
initial conditions and follow from the spectrum of an 
Hermitian Hamiltonian operator obtained by a similarity 
transformation of the Fokker-Planck operator governing 
the evolution of the probability density In general, 
the long-time relaxation behavior is dominated by the 
lowest, i.e. the smallest in its absolute value, nonzero 
eigenvalue of the Hamiltonian, defining the spectral gap. 
As we proceed to show, this is not always the case. 
The relaxation pattern to an equilibrium distribution 
may contain additional terms, with rates that depend 
on the initial distribution and are not present in the 
Hamiltonian spectrum. The smallest non-spectral rate 
can be smaller than the smallest spectral relaxation rate 
and may thus dominate the relaxation behavior over 
the whole observable range. In this case the knowledge 
of the spectrum of the Hermitian counterpart of the 
Fokker-Planck operator does not contribute much to the 
time-dependent solution of the Fokker-Planck equation. 
In this Letter we consider the exactly solvable examples 
of Ornstein-Uhlenbeck processes (OUP) describing the 
coordinate of an overdamped harmonic oscillator driven 
by a Gaussian white noise, and a generalized Levy 



OUP driven by symmetric, white Levy noise. Even 
for the standard OUP the non-spectral relaxation of 
broad initial distributions has so far surprisingly been 
overseen in the textbooks and the applied mathematical 
literature. It was never considered for the Levy OUP, 
for which the similarity transformation to a Hermitian 
operator was not known. The effect of non-spectral 
relaxation found in this letter is of different nature 
than the non-spectral relaxation in the presence of 
multiplicative noise described in 

The time dependent probability density p{x, t) for 
a Gaussian diffusion process in a one-dimensional 
potential U(x) solves a Fokker-Planck equation [1[ of 
the form 

d d d 2 

— p{x,t) = — [U'(x)p(x,t)] + -—£p(x,t) = L v p{x,t), 

(1) 

where the time has been scaled to units of the inverse 
diffusion constant. The Fokker-Planck operator L 1 ^ may 
be indexed by the potential U and the order fi = 2 of 
the derivative in the noise term. The time independent, 
stationary solution is given as p s t(x) — ^e~ u ( x ^ where 
Z is determined by normalization. The time evolution of 
the transformed function 

^(x,t)=p(x,t)/^/p^(x)^ei u ^p(x,t) (2) 

is given by a Hermitian Operator H 



V ^ ~ ft? 



iP(x) (3) 
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with V(x) 



[\U'{xf 



U"(x)} [1]. Equation © tic function p(k,t) = E p [e ikx ]. There it has the form 



defines a similarity transformation between Ly and H 
and hence a transformation ip — Sp from the space of the 
solutions of the Fokker-Planck equation to the space of 
the solutions of the Schroedinger-like equation ©. The 
eigenvalues — A of the Hamiltonian H are real valued and 
the corresponding eigenfunctions ip\(x) form a basis in 
the space of square integrable functions. A distribution 
p(x, t) that can be expanded into the transformed 
eigenfunctions <p\(x) = il)\{x)y/p s t(x) will relax at rates 
that are typically given by the eigenvalue spectrum of 
H. We call this a spectral relaxation pattern. However, 
from Eq.Q follows that only those distributions p(x,t) 
transform into a square integrable function ip{x, t) which 
decay sufficiently faster at infinity than 1/ \/pst(x) 
grows. If U{x) goes faster to infinity than logarithmi- 
cally, p(x, t) must decay exponentially. In this case, the 
existence of all moments is a necessary and, indeed, 
sufficient condition for spectral relaxation. Other, fully 
legitimate probability density functions, for instance 
a Cauchy distribution p(x 7 t ) — (l/n)(x 2 + as 
initial distribution for the Fokker-Planck equation, 
are mapped to functions that grow rapidly at infinity 
and cannot be expanded into the square integrable 
eigenfunctions of the Hamiltonian. The relaxation for 
such initial conditions does not have to be spectral. 
Since one is usually interested in the Green's function of 
the system, which is the conditional probability density 
p(x,t + r\xo,t), square integrability of the transformed 
initial condition is always assumed and other cases 
have never been considered. For the OUP with a linear 
restoring force, given by a mobility coefficient v 1 the 
potentials are U{x) = \vx 2 and V(x) = ju 2 x 2 — \v. 
The spectral relaxation rates are given by the energy 
eigenvalues — A„ = nu, n G N of the quantum harmonic 
oscillator with ground state energy zero. 

Let us proceed to show that the fractional Fokker- 
Planck operator of the Levy OUP [1, 0| can also be 
mapped to the Hamiltonian of the quantum harmonic 
oscillator. The Fokker-Planck equation for the probabil- 
ity distribution of a Levy Flight in a harmonic potential 
reads 



t)p 
di 



dx 



[uxp{x, t)]+ A^ 2 p(x, t) = LZp(x, t) , (4) 



with < /i < 2. The fractional Laplacian is defined by its 
action in Fourier space: A fJ/ ^ 2 p(x) — > — \k\ fJ, p(k), where it 
is diagonal. We write Z/£J for the corresponding fractional 
Fokker-Planck operator depending on the noise param- 
eter fi and the mobility v. Eq.([T]) with U(x) = \vx 2 
is a special case of Eq.(|4]) with /i = 2. The relaxation 
rates are eigenvalues of the operator L£J. In Fourier 
space Eq.(j3]) is an evolution equation for the characteris- 



^p(k, t) = -vk^p{k, t) - \k\»p(k, t) (5) 

By simply rescaling the argument with a di- 
agonal transformation T a with integral kernel 
T a {n, k) = <5(|K| 1 / Q sign(K) - fc) 



T a p 



(M) 



T a (n, t)dk = p{\n\ <* sign(/-c), t) 



(6) 

and using the chain rule kdk[T a p] = ctnd K \T a p] we find 
that, with a = /i/2, the transformed functions follow a 
non-fractional Fokker-Planck equation 



d ~ fi d - 

— TVp = — V — K-— \ThlV - 



Ok 2 



[f^p]^L 2 [f^p]. (7) 



For any a > the transformation is defined everywhere, 
it preserves the value p(k = 0) = [T a p\{n = 0), i.e., the 
normalization in coordinate space, and it has T^ 1 — T±/ a 
as the inverse. Indeed, it defines the similarity trans- 
formation T.. loLf.T^L = L 2 ,„ of the fractional Fokker- 

Planck operator to that of the non-fractional OUP 
L 2 ^^ with the coefficient of restoring force depending on 
the noise parameter /i. In coordinate space, the transfor- 
mation T a is an integral transform with the kernel 



T a (x,x) 



1 

2^ 



£ iKx-i\K\ a sign(K)x^ K 



(8) 



Thus, applying transforms S and T M / 2 m sequence, 
one can transform the fractional operator to the 
Hamiltonian of the quantum harmonic oscillator with 
the harmonic eigenvalue spectrum — A n = nv^r. 



The spectrum of a Hermitian operator and its eigenfunc- 
tions are determined by the properties of the Hilbert 
space it is operating on. In the case of square integrable 
functions, there are selection rules that constrict the pos- 
sible eigenvalues, and the corresponding eigenfunctions 
form a complete basis. Given the similarity transforma- 
tion between the Fokker-Planck operators for the one 
dimensional Gaussian diffusion process in any potential 
or the generalized Levy OUP and a quantum mechanical 
Hamiltonian, it is tempting to use the same selection 
rules in order to also resolve the identity in the space 
L 1 of integrable solutions of the initial Fokker-Planck 
equation. But we understand that this is only possible 
in a subspace of L 1 . However, both, the solution of the 
eigenvalue problem, and the complete time dependent 
solution of the fractional Fokker-Planck equation for 
the Levy OUP, can be found analytically. We therefore 
use this analytically tractable example as a showcase to 
explain non-spectral relaxation, which could otherwise 
not be explained by Hermitian spectral theory. 
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Since we only consider real valued functions in co- 
ordinate space, we can restrict the eigenfunctions (fi\(k) 
of the Fokker-Planck operator in Fourier space to 
those for which (f\(— k) = <p\(k)* holds. The eigenvlaue 
problem L^ipx — Xip\ is solved via seperation of variables 
by any A 6 R, A < and 



<p\(k) = [a x + z&Asign(fc)] |fc| « e 



-Jr-IAI- 



(9) 



with a\,b\ £ HL The restriction of A to non-positive real 
values follows from the existence of a stationary limit 
lim t _ s . 00 p{k 1 1) = ifa(k) and the transformation of to 
a Hermitian operator. A nonzero coefficient a\ means 
that the eigenfunction in coordinate space has a nonzero 
even part and a nonzero b\ contributes to the odd part 
of (p\(x). For the symmetric Levy flight in a symmetric 
potential, the stationary solution must be even, and we 
find Pst(k) = <po(k) = exp ( — I ^ I p ) an d thus 

ipx(k) = [a x + ib x sign(k)} \k\»p st (k). (10) 

The characteristic function p(k,t + r) at time t + r 
is the unique solution of the fractional Fokker-Planck 
equation ([5]) with given initial characteristic function 
p{k,t) — po(k) at time t. It is found by the method 
of characteristics and yields 



p{k,t + r) 



Po(ke~ 



Pst{ke~ 



:Pst(k). 



(11) 



Comparing Eqs. flTU)) and (fTTj) we see, that p(k, t + r) has 
an expansion into eigenfunctions of L£ if the ratio po /p s t 
can be expanded as 



Po(ke-^) 
Pst(ke- vr ) 



^2[a x + ib x sign(k)]\k\ « . 



At 



(12) 



In more general cases, the sum may be replaced by an 
integral with respect to an appropriate measure over the 
non-positive real numbers. Note, that both, the initial 
distribution and the stationary distribution determine 
the relaxation rates to the equilibrium. Two-point cor- 
relation functions at equilibrium require the conditional 
probability distribution p{x,t + T\x ,t) which is the 
solution of the Fokker-Planck equation with a delta 
distribution po{x) — S(x — xq) at an initial time t. The 
asymptotic relaxation rate of these correlation functions 
can be different from the relaxation rates of probability 
densities with other initial distributions. 

In this letter we do not study the conditions under 
which such an expansion exists. Instead, we show an 
example for which the expansion into an absolutely 
convergent series is known, and where the contributing 
eigenvalues A are not identical with the harmonic 
spectrum A„ = -ni/|, n G N. Let us consider a 
Levy stable distribution of index a centered around 



a point xq, which has the characteristic function 
Po{k) = exp (ikxo — <r\k\ a ). The fraction (fT2j) has the 
absolutely convergent series expansion 



Pst{ke- VT ) 



— ^ ^ C-hnn |&| u & 



l,m : n— 



with 

Chnn = 77— Tl (^oSign(fc))' (-Cr)™ ( — 

L\m\n\ 
and relaxation rates 



Azmn = —v(l + ma + nfi), l,m,n G N. 



(13) 



(14) 



(15) 



Note, that odd eigenfunctions occur in the expansion 
only for odd I and asymmetric initial conditions xq =/= 0, 
and the smallest eigenvalue of an odd eigenfunction is 
simply given by the deterministic exponential relaxation 
of the mean to its stationary value zero at rate v, 
independently of the noise parameter (i. Non-spectral 
relaxation rates are observed, whenever the initial 
distribution does not belong to the domain of attraction 
of the stationary distribution (as a stable law), i.e. for 
a =/= /i, or /i < 2 and xo ^ 0. A delta distribution at 
the origin, i.e. xq — and a — 0, can be expanded into 
the even eigenfunctions corresponding to the harmonic 
eigenvalues A2« = —nv[i, because this special case be- 
longs to the domain of attraction of all stable laws. The 
corresponding expansion was found in which could 
mistakenly be interpreted as a hint that the complete 
eigenvalue spectrum of the Fokker-Planck operator for 
the Levy OUP is, in fact, harmonic. 

Given the time dependent solution (fTT|) of the Fokker- 
Planck equation, in order to demonstrate non-spectral 
relaxation, one can look at relaxation of the expected 
values of appropriate observables to their stationary 
values. Instead, here we use the square norm of the 
difference as distances 



A 2 + 



2tt 



{p + {x, t) - p st (x)) 2 dx 



{p + {k lT )-p st {k)fdk, 



(16) 



/OO \ f°° 

p~(x,T) 2 dx = — [p'{k,r)\ 2 dk, 
-oo ^ J — OO 



between the even and the odd parts of the time de- 
pendent distribution p(x,r) — p + (x,r) +p~{x,r) and 
the corresponding parts of the stationary distribution 
Pst(x) = ptt( x )> which has no odd part. The relaxation 
rates of these L 2 distances assume twice the value of the 
eigenvalues in the expansion fj 12|) of po/p s t- As a par- 
ticularly striking example, in Fig[TJ we plot the L 2 dis- 
tances for the case of the standard OUP (/i = 2) with a 
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Figure 1: (color online) Semilogarithmic plot of the L 2 dis- 
tances A+ and A 2 _ between, respectively, the even and odd 
parts of the time-dependent and the stationary probability 
density, in case of the OUP with Gaussian white noise (y = 1, 
fi = 2) and shifted, tempered a-stable initial distribution 
fEq, (|17[l . a = 1, a = 2/3, xo = 1). The asymptotic relaxion 
rate of the square distance of the odd part (thin solid, blue 
line) is — 2Ai = 2v — 2, independent of the cut-off parameter 
7. The square distance of the even part displays a cross-over 
from slow, non-spectral decay at a rate 2a = 4/3 to spectral 
relaxation at rate — 2A2 = 2fi = 4. The transient is longer 
for smaller values of 7, i.e. broader distributions. Here we 
have plotted A+ for 7 = 0, 10~ 6 ,10~ 4 and 10~ 2 (bold, solid 
lines). Using a Levy stable distribution as initial condition, 
i.e. 7 = 0, non-spectral relaxation of A+ is observed at all 
times. The dashed lines are exponential functions | exp(Ar) 
with A = —4, —2 and —4/3, drawn for comparison. The L 2 
distances where calculated according to Eqs. (|ll[) . (|16[) and 
(|17|) by numerical quadrature in Fourier space. 



is delayed during a transient that depends on 7, and 
may be considerably slower than Av, depending on the 
index a of the initial distribution. In fact, the transient 
decay rate is given by twice the smallest eigenvalue 
A = —a used in the expansion of the a-stable law 
approximated by the initial distribution. For very small 
7, before entering the asymptotic regime the even L 2 
distance A3, can become so small, that in experiments 
the cross over may not be observable at all. While 
non-spectral relaxation is a transient phenomenon in the 
standard OUP for broad a-stable initial distributions, 
the eigenvalues A/„ = — v(l + n/i), !,neN, used in the 
expansion of the conditional probability density, with a 
delta distribution as initial condition, are non-spectral 
for any /1 ^ 2, i.e. for a generalized Levy OUP. 

In conclusion, we have shown that the spectrum of 
the Hermitian counterpart of a Fokker-Planck operator 
corresponding to a Gaussian diffusion process in a 
potential only determines the time evolution of suffi- 
ciently localized probability densities. Even in physically 
realistic situations, when all moments exist, broad initial 
distributions may, during a possibly long transient, relax 
slower than expected from the Hermitian eigenvalue 
spectrum. We have also shown that the fractional 
Fokker-Planck operator for a Levy flight in a harmonic 
potential is related to the Hamiltonian of a quantum 
harmonic oscillator. The finding of such a similarity 
transformation gives new insights into the structure of 
the solution space. However, the transition probabilities 
of the generalized Levy OUP, which describe the equi- 
librium statistics for single trajectories of the Markov 
process, cannot be expanded into the transformed 
eigenfunctions of that Hamiltonian. While Hermitian 
operator spectral theory is a powerful tool to analyze a 
system, it is important to understand the limitations in 
its theoretical and experimental applications. 



smoothly tempered, Levy stable initial distribution [5|, |6[ 
shifted to a point xq 0. The characteristic function for 
the jump size distribution of the truncated Levy flight 
with exponential cutoff was found in Q and is given by 

( -i Re[( 7 + ^|) Q ]-7 Q \ 
poO)=exp [tkxo-a — (17) 

V C0S (<* 2 ) / 

Since for 7 ^ all derivatives at k = 0, and hence, 
all moments exist, we expect spectral relaxation at the 
asymptotic rates — 2Ai = 2v for the odd L 2 distance A?_, 
and — 2X2 = Av for the even L 2 distance A^_. In FigUwe 
observe that At relaxes at the spectral rate — 2Ai = 2v 
independently of 7. On the other hand, the decay of A^, 
which is expected to be faster than its odd counterpart, 
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